Definition and relevance of nonequilibrium intensive thermodynamic parameters 
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We show that intensive thermodynamic parameters associated to additive conserved quantities can 
be naturally defined from a statistical approach in far-from-equilibrium steady-state systems, under 
few assumptions, and without any detailed balance requirement. It may apply, e.g., to dissipative 
systems like granular gases where volume or mass is still conserved, or to systems with periodic 
boundary conditions where fluxes of conserved quantities are present. We emphasize the usefulness 
of this concept to characterize the coexistence of different nonequilibrium phases, and discuss the 
influence of the contact between two different systems, in relation with measurement issues. 
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One of the essential features of intensive thermody- 
namic parameters (ITPs), like temperature, pressure or 
chemical potential, is that these parameters take equal 
values in two different systems that can exchange the con- 
jugated extensive quantity (energy, volume, or particles) 
once equilibrium is reached. This is a strong property 
since it holds even if the systems put into contact have 
very different microscopic dynamics. This equilibration 
between different systems makes ITPs a powerful concept 
to describe the influence of the environment or phase co- 
existence, and is essential to allow, say, the temperature 
of a system to be measured with a thermometer. 

Several approaches have been developed in nonequilib- 
rium statistical physics to generalize the notion of ITP, 
usually focusing on temperature. For aging systems, 
an effective temperature may be defined from the long- 
time slope of the fluctuation-dissipation relation (FDR) 
IH 0- For systems weakly driven into a stationary 
nonequilibrium state, ITPs are defined locally as in equi- 
librium Q. On the other side, phenomenological exten- 
sions to far-from-equilibrium regimes seem to lack mi- 
croscopic foundations. For such systems, statistical tem- 
peratures have been proposed in specific models, from 
FDR [|, maximum entropyprinciple 0], or generalized 
microcanonical approaches |8|]. Yet, most of these stud- 
ies do not discuss the relevance and generality of the so- 
defined temperature. 

In this Letter, we propose a general definition of ITPs 
conjugated to additive conserved quantities in nonequi- 
librium systems, valid under few assumptions. This def- 
inition, which does not require any detailed balance re- 
lation to be fulfilled, is illustrated on a stochastic model 
in which a dissipated energy is coupled to a globally con- 
served mass. The usefulness of nonequilibrium ITPs to 
characterize the coexistence of nonequilibrium phases is 
outlined, on the example of a two-species mass transport 
model. Finally, we discuss the problem of the contact 
between two nonequilibrium systems with different mi- 
croscopic dynamics, as well as the relation with the large 
deviation functional approach ||. 



Framework and hypotheses.- Starting with a general 
formulation of the problem, we consider the steady state 
of a macroscopic nonequilibrium system for which some 
additive quantities Qk, k — 1, ...,£, are globally con- 
served by the dynamics. On general grounds, the station- 
ary probability to be in a microstate a, with probability 
weight f a , is of the form 



P a = 



fa 



Z{Qi, Qe 



(i) 



fc=i 



where the product of delta functions accounts for the 
conservation laws, and Z{Q\,...,Qi) is a normaliza- 
tion factor, named hereafter a "partition function" . 
Splitting the system into two subsystems S\ and S2 , 
so that a = {0:1,0:2}, we introduce the probabil- 
ity ^(Qn, . . . , Qn\Qi, . . . , Qi) that conserved quantities 
have values Q k \ in Si, given their total values Q k 

e 

*({Qki}\{Qk}) = n we - (2) 

ai,a 2 fe=l 

The key assumption in the following derivation is that 
this conditional probability satisfies an asymptotic fac- 
torization property, namely 

ln*({Qki}|{Qk}) = ^Z 1 ({Q kl })+lnZ 2 ({Q k -Q kl }) 
- \nZ{{Q k }) + e N {{Q kl },{Q k }) (3) 

where Z v (\Q kv \) refers to subsystem S v , when isolated, 
and €N{{Qki},{Qk}) becomes negligible with respect to 
In Z v in the large N limit, N being the number of degrees 
of freedom. Although not obvious in nonequilibrium sys- 
tems where long-range correlations might develop, this 
factorization property actually holds in some rather large 
classes of models. It is satisfied, for instance, in lattice 
models like the ZRP [T^j, the ARAP or more gen- 
eral mass transport models 0, ^| , where the probabil- 
ity weight f a factorizes as a product of one-site weights 
fa = T\i=i fi,an with a = {ai}. In this case, ejv = in 
Eq. ©. Besides, this asymptotic factorization property 
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also holds in one-dimensional models with conservation 
laws, for which the stationary probability is given by a 
matrix product ansatz with finite-dimensional matrices: 



f a = TrfJ i=1 M ai , where M ai is a matrix associated to 
the local state «i (for nonperiodic systems, one rather 
has f a = {W\l\f =1 M ai \V), with vectors (W\ and |V) 
determined by the boundary conditions). This approach 
has been used, e.g., in the context of the ASEP model, 
and can lead to finite-dimensional [lil Hif or infinite- 
dimensional matrices fl5t fl6j |. We shall simply sketch 
here the argument leading to Eq. (JHJ), and defer a more 
detailed presentation to a later publication ^3] • For sim- 
plicity, we consider a generic model with local variables qi 
and one conserved quantity Q — 5Z»=i 1i- The stationary 
distribution P({qi}) reads 

P({Qi}) = ^n M (*)) 6 (l> - Q ) ' W 

where M(q) is a finite-dimensional matrix. Introducing 
the matrix R(Q) such that Z(Q) = Tri?(Q), one has 

1 



*(Qi|Q) 



Z{Q) 



Tr[Ri(Qi)R 2 {Q - Qi)] (5) 



Denoting by R(s) and M(s) the Laplace transforms of 
R(Q) and M(q), one has i?(s) = M(s) N . If M(s) is in- 
vertible, M(s) = exp(S(s)), so that R(s) = exp(NB(s)). 
B(s) may be decomposed as B(s) = D(s) + L(s) where 
D(s) is diagonalizable, L(s) is nilpotent (i.e. L(s) p ^ — 
0), and D(s) commutes with L(s). It follows that R(s) = 
exp(ND(s)) exp(iVX(s)), where exp(iVL(s)) is a polyno- 
mial in NL(s): the dominant contribution to R(s) is pro- 
portional to N^' 1 exp(JWi(s)), with £i(s) the largest 
eigenvalue of D(s). If p(s) is bounded when s — > 0, one 
can show performing an inverse Laplace transform that 
the asymptotic factorization property © holds 01 ■ 

For infinite matrix representations, this argument may 
break down as the eigenvalues of D(s) might not be 
bounded. From a physical viewpoint, it is reasonable 
to think that the factorization property J2J is related to 
the presence of a finite correlation length in the system, 
as large subsystems then become essentially independent. 
Both properties are expected to hold for finite matrices 
|l5|. Yet, for infinite matrices, whether these two prop- 
erties are actually related or not remains an open issue. 

Definition of nonequilibrium ITPs.- To define a 
noncquilibrium ITP, we compute, guided by equilibrium 
procedures, the derivative of In ^({Qki}\{Qk}) with re- 
spect to Qki at the most probable value Q* k \- Equating 
this derivative to zero yields 



d\nZ x 



dQ 



kl 



d\nZ 2 



dQ 



1,2 



Qk-C 



(6) 



Thus it is natural to define the ITP of the system as 

dlnZ 



At 



dQ k 



(7) 



since this parameter, when defined within subsystems 
that can exchange the quantity Qk, takes equal values 
in both subsystems once the steady-state is reached, as 
seen from Eq. 10 -see for a more detailed discussion 
on a specific example. In the following, this property is 
loosely called the "equilibration" of subsystems, although 
we deal with noncquilibrium situations. 

The above partition into two subsystems is also useful 
when one subsystem is small, but still macroscopic. The 
effect of the rest of the system, acting as a reservoir of Qk, 
may then be encoded into the parameters Xk ■ Integrating 
the distribution over the degrees of freedom of the 
reservoir leads for the subsystem to 



fa 



exp 



k=l 



(8) 



which may be called a "grand-canonical" distribution. 
Note also that In Z may be interpreted as a nonequilib- 
rium thermodynamic potential, and that, as in the equi- 
librium formalism, the cumulents of Qk are given by the 
successive derivatives of In Z with respect to A^ |l7j • 

Illustration on a dissipative model - Let us illustrate on 
a simple nonequilibrium model the definition J7J of the 
ITP. Wc consider a one-dimensional model with two dif- 
ferent physical quantities, say, an "energy" E — ^TJili £ i 
and a "mass" 

M = EiLi™i- Ener 

gy is flowing through 
the boundaries, and is partially dissipated in the bulk, 
whereas mass is conserved and cannot flow through the 
boundaries. The continuous time dynamics is defined as 
follows. An amount of energy cj may be either injected 
at site i = 1 with rate J(u>), moved from site i to i + 1 
with rate ipi(u>\ei,mi) (if i = N, the energy flows out of 
the system), or dissipated on site i with rate A(w|ei, mi). 
Besides, an amount of mass /i may be transported from 
site i to j = i ± 1 (except for i = 1 or N, where j = 2 or 
N—l respectively) with rate tp2{lA £ ii m i)i an d with equal 
probability for both target sites. To our knowledge, this 
model, which generalizes the cascade model introduced 
in [lsj . was not considered previously in the literature. 
In order to allow for a factorized distribution, the rate 
functions arc chosen as follows 



J(w) = «iHe- b " 



tpi(u\e,m) = Vi(lo) 
(p 2 (n\e,m) = v 2 (fi) 



g(e,m) 
g{e,m- p) 
g(e, m) 



(9) 
(10) 

(11) 



A(«|e,m) = (e^-l)^)^^ (12) 

g{e,m) 

with positive functions vi(u), 1)2(1-1), g(e,m), and two 
positive parameters a and b. Clearly, due to the presence 
of flux and dissipation, the above dynamics cannot satisfy 
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detailed balance. The steady-state distribution reads 



Uf=ifi( £ j, m j) 
Z(M) 




M\ (13) 



is Z e (Mi c , M 2c ) — g(Mi c , M 2c ), and the chemical poten- 
tials A£ = d In Z c j dMkc read 



M 2c 



(15) 



with fj(sj,m,j) = g(ej,rrij) exp[— (aj/N + b)ej] on site 
j. Let us choose for the function g(e, to) the simple form 
g{e, to) = exp(— urns' 7 ) with er > 1 and k > 0. The parti- 
tion function can be computed in the large N limit, and 
one finds Z(M) = AM (i-V^-i with A independent 
of M. This gives for the ITP A = (a - l)/(ap) with 
p = M/N, showing that A depends, through a, on the 
coupling to the nonconserved energy. 

Relevance of ITPs for phase coexistence - Let us now 
discuss the case where two systems with the same micro- 
scopic dynamics, but different macroscopic states, can 
exchange some globally conserved quantities. A well- 
studied example of such a situation is the condensation 
transition observed in ZRP fiolll^. Eol Elf , or more gen- 
eral mass transport models 13]. This transition occurs 
when the overall density exceeds a critical value p cnt , 
so that a finite fraction of the total mass condense onto 
a given site. In standard approaches 0], the critical 
density is obtained by looking for the convergence radius 
of a formal grand-canonical partition function [2^]. De- 
tails about the condensed phase need to be studied in the 
canonical ensemble where the total mass is fixed 0]. 

In the following, we interpret the condensation as a 
phase coexistence and show that enforcing the equality 
of the ITP defined in Eq. (JJJ -here, a chemical potential- 
in the two phases leads to a natural quantitative descrip- 
tion of the condensation. To this aim, we compute sep- 
arately the chemical potential of the condensate, taken 
as isolated, and the equation of state of the fluid phase. 
Let us illustrate this point on a concrete example. We 
consider a variant of the model defined by Eqs. JSJ to 
(|12fl . defined on a periodic lattice, where Ei is now a con- 
served quantity (that is, J(uj) = and A(cj|e, to) = 0), 
so that we shall replace the notations Si and rrii with mu 
and mu respectively (both quantities are now "masses" ) . 
The continuous time dynamics proceeds by transferring 
a mass p\ or p, 2 from a random site i to i + 1, with rates 
(£>i(/X;l|toh, TO2i) or ip 2 (p 2 \mu, TO2i) given in Eqs. (JJOJ and 
(|TTJl respectively. In dimension d > 1, the mass p\ or p 2 
may also be transferred in any direction transverse to the 
flux, with equal probability. The steady-state distribu- 
tion in arbitrary dimension d is given by 



P({m ki }) 



lL=lff<>l*> TO 2 4 ) -Q s 



Z{M U M 2 



k=l 



n (14) 

Guided by the two-species ZRP studied in [21J , we choose 
g(mi,m 2 ) = exp(— k m 2 rri~~ v ), with v > 0. Assuming 
that a condensate containing large amounts M\ c and M 2c 
forms, the partition function of the isolated condensate 



Thus in the thermodynamic limit, X 2 = 0. Enforcing the 
equality of the chemical potentials in both phases A^ = 
A£ = Afe , the densities p k in the fluid phase are computed 
in the grand-canonical ensemble, yielding p{ = (^+1) /Ai 
and p{ = T(l + 2i/)/[T(l + v)k\\\. Hence Ai > 0, and 



M lc 



and p[ 



M. 
f _ 



1/(^+1) 



2 c 
Pi- 



The condensate M\ c is subextensive, 



It follows that Ai, and thus p 2 , are fixed 



by the overall density p\ . A condensate forms if p 2 > p 2 , 
meaning that p^ is the critical density p 2 rtt - One has 
M 2c = N{p 2 - pf lt ), and Eq. O for Af gives 



(16) 



Accordingly, the framework of ITPs provides a simple 
description of the condensation in terms of the (out-of- 
equilibrium) "equilibration" of two coexisting phases. 

Contact of systems with different dynamics - In this 
part, we consider the issue of the equalization of ITPs 
when two systems with different microscopic dynamics 
are put into contact. This is precisely the type of prob- 
lem encountered when one wishes to define a "thermome- 
ter" , which in the present more general nonequilibrium 
context, we shall call "ITP-meter". For such an "ITP- 
meter", one requires two essential properties. First, its 
ITP must equalize with that of the system over which 
the measure is performed, without perturbing this sys- 
tem. Second, one needs to know the equation of state of 
the "ITP-meter" to relate its ITP to a directly measur- 
able quantity. This is a highly non trivial problem, which 
depends on the microscopic dynamics of each system, but 
also on the dynamics at the contact. Let us imagine that 
two different systems S\ and £2 , that separately conserve 
the same physical quantity Q = J^jfft, are put into con- 
tact. The dynamics at the contact indeed imposes the 
distribution $(Qi|Q) for the random partition of Q into 
Qi and Q 2 — Q — Q± over the two systems. The probabil- 
ity distribution of the whole system then reads, assuming 
for simplicity that the probability weights of each system 
strictly factorize (the index i v refers to S u ) 



ME* 



Zi(Qi)Z 2 (Q - Q x ) 

(17) 

with Qi = SiieSi^i- The ITPs of the two systems 
equalize only if $(Qi|Q) = Z x {Qi)Zi(Q - Qx)/Z{Q). 
This relation is indeed satisfied in equilibrium due to the 
detailed balance relation, if the Hamiltonian is additive 
|17| . Out of equilibrium, the above factorization and the 
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equalization of the ITPs does not hold in general, but are 
recovered inside some classes of nonequilibrium systems. 

To illustrate this point, let us consider two single 
conserved mass transport models Si and 5*2 with site- 
independent transport rates ^1.2 (/i|™) = w(m)/i,2(w ~ 
A*)/ fi,2(m), with different fi(m) and /2(m), but the same 
v(p). The dynamics at the contact is defined as follows: 
a mass leaving a boundary site may be transferred to 
one of the neighboring sites or, with the same probabil- 
ity, to the other system. The global probability weights 
factorize (as v(/J,) is site-independent 0, 0|), and the 
ITPs of the two systems equalize. The key point is that 
the equalization of the ITPs occurs for two different mi- 
croscopic dynamics (fi(m) ^ j ^(m)), but belonging to a 
given class (same v(p)). Indeed, if Vi(fJ,) ^ w 2 (^), such an 
equalization does not hold. This result suggests the exis- 
tence of classes of systems which mutually "equilibrate" , 
so that any member of a class, with known equation of 
state, may be used as an "ITP-meter" for the other mem- 
bers, under suitable size ratios. 

Discussion- We have proposed a general definition 
of ITPs associated to additive conserved quantities in 
nonequilibrium steady states, under the asymptotic fac- 
torization assumption Q. This property holds in par- 
ticular if the stationary distribution is factorized or, for 
one-dimensional systems, if it can be computed from a 
matrix product ansatz with finite-size matrices. Its ap- 
plicability to more general situations remains an impor- 
tant open issue. Indeed, it would be essential to have a 
criterion allowing to determine experimentally or numer- 
ically whether the factorization property holds, without 
having to know the stationary probability distribution 
explicitely. Note that in equilibrium, the same difficulty 
arises in principle, but it reduces there to the "additivity" 
of the hamiltonian (short-range interactions). 

Besides, we note that if ^(Qi\Q) may be written as 
*(Qi|Q) = exp[-NG(pi,p)], with pi = Qi/N x and 
p = Q/N, Eq. J2J corresponds to an additivity prop- 
erty for G{pi,p), which somehow resembles the additiv- 
ity principle for the large deviation functional J-({p(x)}) 
of the density profile p{x) in the open ASEP model 
As the physically relevant states minimize G(pi, p) with 
respect to pi in the thermodynamic limit, G{pi, p) ap- 
pears as a "free energy" function |2jj- In analogy with 
the ASEP where !F({p(x)}) is a nonlocal functional of 
p{x), G(pi, p) is a nonlocal function as it depends on the 
global density p, due to the conservation law. Still, the 
physical origin of this nonlocality may be different in the 
two problems. Also, the quantity Q is conserved here, 
whereas particles are exchanged with boundary reservoirs 
in the ASEP. 

From Eq. ©, one sees that the equality of ITPs be- 
tween subsystems is related to the minimization of a non- 
local "free energy" function. Note however that this is 
already true in equilibrium systems. The specificity of 
nonequilibrium dynamics shows up mostly when the sys- 



tems into contact have different microscopic dynamics. 
In this case, the ITPs of the two systems may not equal- 
ize, due to the essential role played by the dynamics at 
the contact -a role hidden at equilibrium due to detailed 
balance. Yet, there may exist classes of nonequilibrium 
systems that mutually "equilibrate" . 
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By "formal" , we mean that it is defined by analogy with 
equilibrium as a superposition of different canonical par- 
titions functions, without being clearly associated with a 
well-defined physical ensemble -see, e.g., [T^ . 
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tems, one may also construct in this way a large deviation 
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